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Abstract 

We propose a method to measure the Witten index using lattice simulation. A 
requirement for the lattice model is that it has at least one exact supersymmetry at 
finite lattice spacing. We prove the validity of the method in case of the supersymmetric 
quantum mechanics, where the index is well known. 

1 Introduction 

Recent developments on lattice formulations of supersymmetric theorie^ make it possible 
to perform numerical simulations to study non-perturbative aspects of the theories. The 
most important feature in these developments is that the action has (at least) one exact 
fermionic symmetry at finite lattice spacing. This was pioneered in [2] for gauge theories 
and in [3] for non-gauge theories. It turned out later that this exact symmetry is the 
scalar part of supersymmetry in terms of the topological twist. The symmetry is strong 
enough to guarantee an automatic restoration of the remaining part of the supersymmetry 
without any fine tuning for the 1- or 2-dimensional case. This restoration is manifestly 
confirmed numerically in 2-dimensional M = (2, 2) super Yang-Mills j4] using a model by 
F. Sugino [5], which also has one exact supersymmetry at finite lattice spacing. A more 
ambitious approach which tries to keep all of the supersymmetry was also proposed by 
D'Adda et.al. in ^ and its Hopf algebraic structure is studied in \7\. A new proposal along 
this line is found in [5]. However, without introducing non-locality and/or extra structure 
such as noncommutativity, only a part of supersymmetry out of full supersymmetry can 
be realized on the lattice; Especially, the equivalence of [6] to other approaches is well 
studied in [5] . Not only simple spacetime lattice regularization but also regularization in the 
momentum space [TU] sheds some light on how to regularize supersymmetric theories (see 
also [TT1[T2]). A different approach based on the large- limit is proposed in [T3]. Further 
developments in the formulation on the lattice are found in [21 USl UHl HZ] • 



*kaiiajnori@to . inf n. it 

^For a recent review, see Ref. [Ij. 
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One of the most interesting non-perturbative aspects of supersymmetry is the sponta- 
neous breaking of the symmetry. Supersymmetry is broken in our current universe and 
the breaking should be the result of a non-perturbative effect, unless it is broken from the 
beginning at the tree level. To study the breaking of supersymmetry using lattice simu- 
lation, a Hamiltonian approach was used in |18) . The present author with H. Suzuki and 
F. Sugino proposed a method for detecting the spontaneous supersymmetry breaking using 
lattice simulation in [12] , where they measured the vacuum energy. They also pointed out a 
possible relation between spontaneous supersymmetry breaking and the existence of a sign 
problem in the simulation, which comes from the phase of the Pfaffian of the fermion bilinear 
operator. 

Another line of application with large- iV limit is found in the gauge/gravity duality. In 
particular, extensive studies of 1-dimensional model for this purpose are found in Refs. j20[ 
W] . The large- iV property of 2-dimensional A/" = (2,2) system has been studied with the 
same motivation [22| . although this system has smaller number of supersymmetries. 

In this paper, we propose another method of detecting the spontaneous breaking of 
supersymmetry. We measure the Witten index [23]. If the index is not vanishing, there 
exists at least one supersymmetric vacuum and supersymmetry is not broken. If the index is 
vanishing, supersymmetry may or may not be broken. Since the Witten index is a partition 
function with periodic boundary conditions in the time direction |24[ 125] . it is of crucial 
importance to determine the normalization factor in simulation. Usually what we obtain in 
the simulation are expectation values normalized by the partition function, but what we need 
here is the normalization itself^ Our idea to determine the factor consists of two ingredients. 
First, we regard the derivation of the path integral formalism from the operator formalism 
as a lattice regularization of the time coordinate. This gives a correct normalization of 
the path integral measure. Second, we measure a suitably chosen quantity, which cancels 
a distribution functional in the path integral. Because of this cancellation, we can 
separate the contribution from the normalization factor. We use supersymmetric quantum 
mechanics (of Af = 2 Wess-Zumino type) [37] as an arena for checking our method. This 
system can be formulated on the lattice in various way^. There is a formulation which 
keeps a nilpotent supercharge Q at finite lattice spacing, with which the action is written in 
a 0-exact form [SHIISIIES- (See also [33].) 

In supersymmetric Yang-Mills quantum mechanics with maximal supersymmetry, which 
would be an interesting application of our method, the index is of crucial importance. This 
system is a candidate for M-theory [31], and to guarantee a suitable limit which gives 
supergravity, the index should be one. Related calculations are found in [35] [3S1 [37] . For 
this system, a numerical treatment in the Fock space is proposed [35] [33 HOI HI] which is 
also useful to obtain the Witten index as well. To the author's best knowledge, however, 
only a little is known about the index from direct calculations of the quantum mechanics so 
far. 

In the next section, we describe our idea in detail. Then we check the validity of the 
measure by analytically calculating the index for supersymmetric quantum mechanics in the 
free case in section|3| In section^, we numerically calculate the index of the supersymmetric 
quantum mechanics and demonstrate that our method in fact reproduces the known index. 

^ An interesting simulation method which allows to determine the normalization and thus the index is 
proposed [17[ on the basis of the Nicolai map. In a 0-dimensional system, Monte Carlo integration was used 
to calculate the partition function in [26 . 

^ Recent proposals with lattice or non-lattice are found in | 10l 1281 1111 1121 18] . See also |29| . 



2 



2 Basic Idea 



The idea for measuring the Witten index is made of two ingredients. One is for the regular- 
ization of the path integral measure and the other is for the regularization of the integraL 

Let us start with the quantum mechanics. The standard way to obtain the path integral 
formulation from the operator formulation is discretizing the temporal direction and then 
inserting a complete set at each of the time slice: 



n 

^ k=l 



dqk 



i){qN- 



-i|e \qN- 



. . n . . ^^-'^ 

where \qi) is a normalized state, H is the Hamiltoniaro of the system, and the time difference 
of the initial state \qini) and the final state \qfin) is T — aN. We regard this a as the 
lattice spacing. Then the standard derivatioi|f| gives the path integral measure for a bosonic 
dimensionless lattice field 6^'^^ as follow: 



. N-l 

/ = n 















/ ° 

' —oo 





For dimensionless fermionic variables, we obtain 



(2.2) 



N-l „ 

vr'D^= n / drkd4>k- 

fc=0 



(2.3) 



In this way, we can determine the natural measure for the path integral. 

In the lattice simulation, what we calculate is not the integration discussed so far. In- 
stead, we calculate ensemble averages which give the ratios of the integrations. We have to 
establish a relation between these averages and the path integral regularized on the lattice. 

In terms of the path integral, the expectation value is 



{A) = 



JV(f>A, 



(2.4) 



where S is the action of the system and we omit fermions for a while. The normalization, 
that is, the denominator is exactly the partition function which we are interested in. Since 
it is made of a ratio and does not depend on the normalization of the path integral, it seems 
impossible to measure the partition function with the correct normalization. However, let 
us consider the following quantity: 



(e 



/P0exp[-1E.M^('^1^^)^] ^ C 



-s ' 



(2.5) 



where fj, is an arbitrary dimensionless parameter. With the measure defined in eq. (j2.2p . we 
obtain 

AT — ^ r . ^ 

d(j)i 



C 



V(j)e 



N-l 

n 

i=0 



27r 



(2.6) 



^ We explicitly denote the operators with hat (") in this section. 

^ For the sake of completeness, we give a brief review of the derivation of path integral formulation in 
Appendix ^ 
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Here, we have renamed the dimensionless lattice field 4)^^ — >■ For any real value of /i, 
we can determine C analytically which gives the overall normalization of the path integral. 
Note that C is given just as an integration defined in (|2.6p . which is independent from the 
distribution of Especially, we never use Gaussian distributed quantities. The left hand 
side of eq. (j2.5l) is an observable in the simulation. Therefore, we can obtain the value of 
the partition function through the following expression: 



Z = / V(j)e 



-s 



C 



p2 0f 



(2.7) 



A comment on the choice of eq. (|2.5p is in order here. We have chosen the Gaussian functional 
g-jIZiM gq^ (j2.5p for simplicity. As one can easily see, it can be replaced with any 

functional of as long as it gives a convergent and calculable value as in eq. (j2.6p . 
Now let us introduce fermions. After integrating out the fermions, we obtain 



-Sf 



V(i)(j[D]e 



-s' 



where 5b and S-p are the bosonic and fermionic part of the action, S' — — InjdetZ?! 
is an effective action with the kernel D of the fermion bilinear, and u[D] is the sign (or 
complex phase) of detulfl The configurations for the ensemble average are generated using 
the effective action S' ignoring the sign factor a[D\. Let us denote the ensemble average 
over these configurations without the sign factor as ( • )o. The sign factor a[D] should be 
reweighted in the measurement afterwards: 



-s' 



f V(j)Aa[D]t 
JV(t>a[D]e-S' 



(2.9) 



To determine the normalization of partition function, we also have to invert the effect of 
the sign factor in addition to the contribution from the effective action e""^ . The analogue 



of eq. (|2.5p becomes {<j[D] 



i) = C/ J 'D(j}a[D\e ^ , where the numerator gives 



again the integration (j2.6p . Now we obtain the expression for the Witten index w: 



w 



c- 



(2.10) 



O.P 



provided all the fields are imposed the periodic boundary conditions as indicated by the 
suffix P. 

Using pseudo fermion ip, we can rewrite eq. (|2.10p as follow. Our definition of the measure 
gives exactly 



'D'ip*'Dip exp 



= j Vtp^^^V tp^^^ a[D] cxp 



-1/2 



where ip = "^('^^^'' + i'P*-^'')- Then, we obtain the expression for the index as 

(^[^pDo.p 



= detD, 
(2.11) 

(2.12) 



'0,P 



Depending on the number of the fermions, the determinant should be replaced by Pfaffian. 
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where 



S"^Sb+Y.^KD^D);//'^,, (2.13) 
C^ = J g-E, = ^-2N_ (2.14) 

It is interesting to see that there is a relation with the sign problem. Eq. (|2.9p imphes 
that the phase quenched average of the sign factor {a[D])o is almost the partition func- 
tion. If this average is close to 0, we cannot obtain reliable expectation values numerically. 
This is the sign problem and in fact it occurs for supersymmetric quantum mechanics with 
spontaneously broken supersymmetry |19) . Note that if the Witten index is vanishing, the 
partition function is vanishing so the expectation value becomes indefinite (or divergent) 
with periodic boundary conditions. Therefore it seems reasonable to regard (not) having 
the sign problem as an indication of the (non-)vanishing of the Witten indexlf] However, no 
justification for this reasoning has been known, because no relation between the index and 
{a[Dp])o^P has been known so far. Even if ((t[Z3p])o,p is vanishing within the small error in 
the simulation, say 0.01 ± 0.02, without knowing the correct normalization, it could mean 
that the index is 0.001 ± 0.002 or 1.0 ± 2.0. Eq. (|2.10p makes a connection between these 
two quantities and we can finally obtain the index from the sign factor. 

One important disadvantage of the method is that it spoils the philosophy of the im- 
portant sampling method due to the factor e^"^ . Therefore, the choice of the value of /j, 
is important. We have to chose it to have as large an overlap as possible between config- 
urations and the operator we measure, namely exp[-|-S" — i ^ - fJ,'^4>i]. The only adjustable 
parameter to obtain larger overlap is fi. 

Finally, it is straightforward to extend our method to higher dimensional cases in prin- 
ciple. The efficiency of the simulation, however, could be far from practical because of the 
overlap problem we have just mentioned above. 

3 Analytical check 

The derivation of the path integral formulation from the operator formulation assumes the 
continuum limit. The lattice artifact may or may not spoil the previous argument. In this 
section, we discuss the effect of the exact supersymmetry on the lattice. We also give an 
explicit calculation of the index for the free case in supersymmetric quantum mechanics. 

Let us suppose that the lattice action has an exact supersymmetry Q, like the one 
proposed in |5] . Then bosonic states and fermionic states make pairs even with finite lattice 
spacings, except for states with Q^Q = 0. At this stage, the situation is exactly the same 
as in the continuum. The difference on the lattice is that we do not have as an exact 
symmetry on the lattice. However, supersymmetric vacuum must be annihilated by Q and 
there is Q-exact Hamiltonian, so the conclusion does not change from the continuum caseH. 
Even with finite lattice spacing, if the index is a non-zero integer supersymmetry is not 
broken. The advantage of exact invariance for the Witten index is also pointed out in |32) . 

A more rigorous argument is the following. Using the fermion number operator F, the 
Witten index is written as 

w = Tr((-l)^e-'3^). (3.1) 

^ The fact that {(t[D])o,p has a close relation to the partition function and thus the Witten index has 
been pointed out in 1191 . See also I42| for a related numerical result. 
^Except for the overall sign. See the discussion below. 
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For simplicity, we assume that the spectrum is discrete and the index does not depend on 
/3 (> 0). For a Q-exact action S = QA with — 0, it is natural to assume that the 
Hamiltonian H in eq. p.ip is also Q-exact: 

H = QJq. (3.2) 

In fact, there is a natural way to define the Q-exact Hamiltonian |19l. 

Under the assumption of the Q-exactness of Hamiltonian, we can show that eq. (j3.1|) 
gives the index even with the finite lattice spacing. The Hamiltonian is complex in general 
in the Q-exact formulation so we have to introduce left- and right-eigen states, which satisfy 

H\Li)=E,\Li), {R,\H ^ {R,\E,, {R,\L,)^S,,. (3.3) 

The Witten index become 



w 



^(i?,|e-^^(-l)^|L.) (3.4) 



in this basis. Since Q commutes with H (note that H ~ {Q, Jq} in the operator notation), 
\Li) and Q\Li) have the same eigen value Ei for H. These two states have opposite statistics 
so they do not contribute to the index, except for the case Q\Li) = 0. Similarly, {Ri\ and 
{Ri\Q have the same eigen value and thus do not contribute except for the case {Ri\Q = 0. 
Therefore, only i such that Q\Li) — and {Ri\Q — contribute to the summation in 
eq. p.4p . These states have Ei = 0, in fact: 

E, = {R\H\Li) = {R,\{QJo + ■hQ)\U) = 0. (3.5) 

Eventually, the index counts the number of the Ei = {) states with a weight factor (— l)-*^, 
which is exactly the same as in the continuum case. 

This result holds even at finite lattice spacing, because we have not assumed anything 
about the continuum limitlfl In particular, the index on the lattice must be an integer. 

We can check the argument above and our measure for the path integral using a trivial 
example, i.e., the free supersymmetric quantum mechanics. The index should be 1 for the 
Q-exact lattice model. 

The action of supersymmetric quantum mechanics with an exact Q symmetry is 

N-l 

S = -QY,i^^k{Fk + d+^k + W'{cj,k)) (3.6) 

k=l3 



N-l 



= E [^('^'^+1 - '^fe)' + l^'i'l'kf + (0fc+i - Mw'i<i)k) 



fc=0 



1 9 

2 



+ M^k+i-i'k) + W"{^k)i>ki^k\- (3.7) 

Here, (f>k is a real scalar, Fk is a bosonic auxiliary field, and tpk is a complex fermion, 
= (t>k+i ~ 0fc, W{(f>k) is a potential and prime (') denotes the derivative w.r.t. (jjk- All 
fields are rescaled to be dimensionless. Q-transformation, which is nilpotent, is: 

g^fc^V-fc, Qi^k=0, (3.8) 

g^fe = Fk - d+(j>k - w'{4>k), QFk = a+Vfc + w"{4>k)i^k- (3.9) 



Of course in order to guarantee that the continuum limit describes the target theory correctly, we have 
to check the restoration of full supersymmetry other than Q. 
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In the following, we assume the free case: 



W{<i>k) 



(3.10) 



where m = amphys is the dimensionless mass. The partition function with the periodic 
boundary conditions is 



Zp = 



n 



dF, 



dip^ dipi I exp 



(detS)"2 detD, 



where 



B 
D 



(3.11) 



(3.12) 
(3.13) 



and A'^-periodicity is assumed for Bij and Dij. In obtaining eq. (13. lip , we made use of 
the replacement F — >■ zF in the integrand to make the integration over the auxiliary field 
F convergent. This gives unity with our measure and the contribution from the auxiliary 
field simply disappears. Alternatively, one can regard that we have moved to the on-shell 
formulation eliminating the auxiliary field, thus the integration should not contain F from 
the beginning which gives the same result. It is straightforward to calculate the determi- 
nants det B and det D. Using the momentum representation, and keeping all the overall 
normalization, we obtain 



det 5 = (4 - 4m + m^)m^ ]J 

N__- 
2 

detD = (-2 + m)m 



(1 - to)4 sin^ — 



k=l 



(1 — m)4sin^ — + 



(3.14) 



(3.15) 



for even TV. Assuming m < 2, which is the case if the lattice spacing is small enough, we 
obtain 

Zp = -l {N: even). (3.16) 

For odd N, we obtain 



det B = 

k=l 

2 

det £> = m ]J 



and thus 



k=l 



Zp 



(1 - to)4 sin^ — + 



(1 — m)4sin^ — + 



{N: odd). 



(3.17) 
(3.18) 

(3.19) 
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Except for the minus sign for even N, we obtain the correct Witten index for the free case 
as expected. 

What is the origin of the extra sign for even N7 It comes from det D. In the derivation 
of the path integral, we have 



{lpk+l\-ip^\ll'l) = V'fc+l'0fc(V'fc+l|'0fc), 



(3.20) 



that is, V'fc+i (='0fc+i) is always combined with ipk- See eq. (|A.23I) in Appendix 1X1 A 
deviation from this combination would give the derivative as in the case in the kinetic 
term. For the action we have used here, however, a natural combination is ipi^ipk- The 
difference seems merely a lattice artifact but it gives a factor nlLo e^'' = (-1)^+^ in the 
determinant, which explains the extra minus sign of det-D. On the other hand, the bosonic 
part does not have such an artifact because 



{qk+i\q\qk) = {qk+i\qk\qk) = {qk+i\qk+i\qk) , 



(3.21) 



so it does not give any extra sign. 

In the case of the naive action, the result is different. The lattice artifact survives and 
gives an extra contribution to the index. Let us start with the following action: 



N-l 
1=0 

Therefore, we have 



_i - 2(j>i) + ^m^(/)^ + tpi{fpi+i - ipi) + rmp^iJi 



(3.22) 



which gives 



det D" 



(detB"a™)2 



B 



'jk 



— 2Sjk — Sj,k+i — ^i+i,fe 



r~\ naive 

^jk 



-2 



Sj+l,k 



^Jk 



a/4 + m2 

Af-l 1 



jv -I 

4(l-m)sin^f 



k=l 



n 

fc=i 



4(1 



4sin2 ^ 



4 sin 



2 far 

TV 



{N: even) 
(N: odd) 



(3.23) 
(3.24) 



(3.25) 



The contribution from the factors in front of the product is =Fl in the continuum limit, where 
TO = oTOphys 0. For even N, the extra minus sign appears again. 

The problem is the inside of the product. Since there is no exact relation between bosonic 
modes and fermionic modes, there is no exact cancellation. Each term has a structure 



I-O(to) - 1-0(1/^^). 
In the continuum limit, we obtain 

det D'laivc 



(det 5" 



N/2 

n(i-^ 



fc=i 



e 2 



(3.26) 



(3.27) 



with a non-zero real number a. It does not give the correct index even in the free case so 
we cannot use the naive action to measure the Witten index. 
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set 


^phys 


set 


\phys 
^2 


^phys 


set 


^phys 


\piiys 
A4 


2a 


0.5 


3a 


4.0 


4.0 


4a 


1.0 


1.0 


2b 


1.0 


3b 


4.0 


16.0 


4b 


4.0 


1.0 


2c 


4.0 


3c 


4.0 


32.0 


4c 


4.0 


4.0 






3d 


2.0 


16.0 









Table 1: Set labels for the free case (n = 2, left), the n = 3 interaction (middle) and the 
n = A interaction (right). We set Lphys — 1- 

From the above argument, if the action has only 0{a^) or higher lattice artifact (even 
without having any part of exact supersymmetry) , presumably we can use it to measure the 
Witten index. For dimension d > 2, the same argument allows only ©(a^"*"^) and higher 
lattice artifact, because the number of the product is ~ N'^. In the case of higher dimensions, 
however, we will need a fine tuning to obtain the supersymmetric continuum limit so the 
argument here is too naive. 

4 Numerical result 

In this section, we check our method using Monte Carlo simulation for supersymmetric 
quantum mechanics on the basis of the lattice formulation (j3.7p . We use the Hybrid Monte 
Carlo algorithm. See [H] for the implementation for this svstenP'l. We fix the physical 
volume iphys = 1- The lattice spacing is a = Lp\^ys/N = 1/N, where N is the number of 
the lattice sites, and the bare coefficients Xi in the potential are scaled as Xf^^^N-'/^. The 
parameters are summarized in table [1] 

4.1 Periodic case 

The very first check is the free case with periodic boundary conditions. There is only one 
bosonic vacuum and the index should be 1 . We list the result in table [5] and plot the index 
versus fi^ in fig. [1] We set fj,^ ~ 0.5, 1.0, 1.5, 2.0, 2.5, 3.0. The index from the simulation is in 
fact 1 within the error (at worst within 3 standard deviations) for odd N and —1 for even 
A'". This result exactly agrees with what we calculated in the previous section. 

As table [5] shows, large |mphys — X2^^^ (set 2c) requires less statistics, because the larger 
mass forces the scalar field to stay around the origin, which makes larger overlap with the 
quantity we measure. Note that large mass corresponds to low temperature measured in 
unit of mphys, since it gives large Lphys'Tiphys- 

Next case is the supersymmetric case with interactions. We use the n = 4 potential 

W{(b) = X2(j)^ + Xicl)\ (4.1) 

In this case, the index is known to be 1 for A2 7^ 0. In table [31 we list the results for the 
index and in fig. [2] we plot the fi^ dependencies. It reproduces the correct index within the 
error (except for /i^ — 0.5). In particular, with a suitable choice of fi^ which gives small 
enough error, we can easily identify the integer. For some values of /i^, 0.5 for example, the 
error is rather large and we cannot determine the integer value. We expect, however, that if 

In addition, we change the time step 5t in the leap-frog precess as follows. In each trajectory we monitor 
the reversibility and if it is broken (it occurs if 5t is not small enough) we restart that trajectory using the 
same initial canonical momentum but smaller Sr. Poor reversibilities are caused by very small eigen values 
of the Dirac operator. 
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Figure 1: The index in the free case. It gives +1 for even N and —1 for odd N as discussed 
in the previous section. 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 



set 4a, N=21 
set 4b, N=21 
set 4c, N=21 



0.5 



1.5 2 

..2 



2.5 



Figure 2: The index for the n — A potential with various values of fjl^. The dashed line 
represents the known result. 



we used larger statistics it should converge to 1. For the strong coupling (set 4c), we found 
that the configuration contained a non-negligible amount of artifact configurations at larger 
lattice spacing. We observed that for some bunch of configurations the artifact surface term 
becomes as much as 10 times larger (or more) than the other part of the action. This is 
the same phenomena reported in [33] in the 2-dimensional system. However, such artifact 
configurations disappear for smaller lattice spacingl^ We observed this phenomena only for 
the n = 4 case. 

In fig. El we plot the dependence on the lattice spacing 1/iV, which shows that the index 
we measured is almost constant within the error against the lattice spacing. 

Odd n in the potential should give spontaneous breaking of supersymmetry. The Witten 
index is known to be 0. We use n = 3 potential: 

W{(b) = A2(/.2 + \3(j)\ (4.2) 

We list the results of the simulation in table S] and plot some of them in fig. [3] and [5l against 

The author thanks Hiroshi Suzuki for pointing out the disappearing. 
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1.04 
1.03 
1.02 
1.01 

1 

0.99 
0.98 
0.97 





set 4b, |x^= 


2.0 1 — ' — 1 











0.02 0.04 0.06 0.08 0.1 0.12 

1/N 



Figure 3: The index for the n = 4 potential with different lattice spacings. 



0.05 


-0.05 
-0.1 

-0.15 
-0.2 

-0.25 
-0.3 



set 3a, N=21 I — h- 

set 3b, N=21 ^-x- 
setSc, N=21 

set 3d, N=21 i B 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 



Figure 4: The index for the n — 3 potential with various values of fi^. The dashed line 
represents the known result. 



IJ, and 1/N, respectively. It gives in fact as expected. 

A comment on the statistics is in order. Since our measurement does not follow the 
philosophy of important sampling, it has an overlap problem in general. Therefore we need 
huge statistics. We use large enough number of independent configurations A^confs to observe 
the behavior of the error ~ -^c^nfs least for some values of fi^ which give small error. See 
also the bottom of fig. [S] for a typical behavior of the error. For some values of /z^ which 
give large error (typically /i^ = 0.5), the statistics are not large enough, since the error did 
not behave as ~ (./Vconfs)"^^^- 

In the measurement, we used configurations every 10 trajectories for n ~ 2, 4 case and 
100 trajectories for n — 3. We used O(10^)-O(10^) configurations which may or may not be 
independent. We estimated the error using the Jackknife method with several bin sizes. For 
n — 2, 4, there were almost no autocorrelations in the measurement so that the bin size we 
adopted is 1 for most sets of the parameter. For n = 3 case, however, the autocorrelation 
becomes longer and the bin size we adopted is 20-250. This is because the sign a[D] does 
not change so frequently in the simulation. The other parts of the observable like S", which 
govern the convergence of error, have almost no auto correlations. Therefore, even if the 
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Figure 5: The index for the n = 3 potential with different lattice spacings. 



set 


N 


num. 




index 


2a 


11 


499000 


2.0 


1.014(9) 


2b 


11 


499000 


2.5 


1.017(15) 


2c 


11 


499000 


1.0 


0.9997(18) 




15 


19900 


1.5 


0.990(7) 




16 


99800 


1.5 


-1.004(6) 




21 


499000 


2.0 


0.989(17) 



Table 2: The free {n—2) case with periodic boundary conditions. We pick up fj,^ which 
minimizes the error from fi^ — 0.5, 1.0, 1.5, 3.0. The num. refers to the independent 
number of configurations after binning. 





1 

e 


et 3d,' ^^=v 


3.0 1 — ^ — 1 




1 







number of independent configurations is smaller than in the n = 4 case, which is governed 

^confs ' 



by cr[D], the number of the configuration is large enough to give the error ~ N^^ 



4.2 Anti-Periodic case 

One may use anti-periodic boundary conditions for fermion in generating configurations. 
Using a reweighting method, we can afterwards obtain the desired index. In this case what 
we measure is 

where we indicated the boundary condition, periodic (P) or anti-periodic (AP) explicitly in 
the sufSx. Note the sign factor a must be calculated from the periodic Dirac operator Dp . 

In fig.[S]and[71 we plot some of the result against various values of ^P. The detailed results 
are given in table [S] It reproduces the desirable result as in the previous periodic case with 
suitably chosen /i^. Note that under the anti-periodic boundary conditions, supersymmetry 
is broken because we impose different conditions for the boson and the fermion. Even 
though, the n = 4 case reproduces the correct index, although the temperature of the 
system is rather low for set 4b thus the boundary condition may not give a large difference 
for the configurations. A rather non-trivial feature is found for the n = 3 case. The average 
of (T[r'p] is not necessarily vanishing but the index is vanishing. For example, (ct[Z?p])o,ap 
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set 


N 


num. 




index 


4a 


11 


199000 


2.0 


1.008(13) 




15 


499000 


2.0 


1.018(19) 




21 


999000 


2.5 


0.88(5) 


4b 


11 


99000 


1.0 


0.999(2) 




15 


49000 


1.5 


1.016(12) 




21 


999000 


2.0 


0.984(12) 


4c 


11 


99000 


2.5 


0.943(10)* 




15 


99000 


1.0 


0.999(6) 




21 


499000 


1.5 


0.989(11) 



Table 3: The n = 4 case with periodic boundary conditions, /x^ is chosen from 0.5,1.0,...,3.0 
to minimize the error. Configurations for the set 4c with N = 11 contain artifact configura- 
tions so the value (*) is not reliable. 



set 


N 


num. 




index 


3a 


11 


1326 


4.5 


-0.008(14) 




15 


4995 


2.0 


0.002(5) 




21 


796 


1.5 


-0.024(23) 


3b 


11 


1237 


4.5 


0.005(4) 




15 


1326 


3.0 


-0.0018(21) 




21 


3118 


4.0 


0.0004(7) 


3c 


11 


4950 


4.5 


-0.001(3) 




15 


2450 


3.5 


0.003(4) 




21 


1237 


3.0 


-0.0009(8) 


3d 


11 


3326 


4.5 


-0.0010(14) 




15 


1243 


3.0 


-0.0006(12) 




21 


660 


3.0 


-0.0005(6) 



Tabic 4: The n = 3 case with periodic boundary condition. is chosen from 0.5, 1.0,. ..,4. 5 
to minimize the error. 
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Figure 6: The index for the n = 4 potential with various values of fi^. The configurations 
are generated with anti-periodic condition. The dashed line represents the known result. 
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Figure 7: The index for the n = 3 potential under various values of /i^. The configurations 
are generated with anti-periodic condition. The dashed line represents the known result. 



set 


TV 


num. 




index 


set 


N 


num. 




index 


4a 


11 


499000 


2.0 


0.999(8) 


3c 


11 


1980 


4.5 


0.002(7) 




15 


999000 


2.0 


0.996(13) 




15 


1960 


4.5 


0.000(3) 




21 


999000 


2.5 


0.98(5) 




21 


6633 


3.5 


-0.0002(6) 


4b 


11 


99000 


1.0 


0.9972(23) 


3d 


11 


1980 


4.5 


0.0007(27) 




15 


99000 


1.5 


0.990(8) 




15 


1237 


2.5 


0.000(1) 




21 


499000 


1.5 


1.012(14) 




21 


1980 


4.5 


-0.0004(4) 



Table 5: The index obtained from configurations with anti-periodic boundary conditions. 
fi^ is chosen from 0.5,1.0,...,3.0(set 4a, 4b) or 0.5,1.0,...,4.5(set 3c, 3d) to minimize the error. 



for set 3d with = 21 is —0.018(7), which is slightly off from 0, while the value of the index 
from fj,"^ = 4.5 is -0.0004(4). The result shows that the reweighting in (|4.3p works perfectly. 

4.3 The less important sampling method 

Since our method contradicts the philosophy of the important sampling, we can prepare 
the configurations using "less important" sampling. We replace the effective action S' in 
generating the configuration by (1 — r)S' , where r is a real number. The ensemble average 
gives (assuming the large enough number of the statistics) 



_ / V<j)A[. 



-(l-r)S' 



In order to obtain the original expectation value, we have to use the reweighting method. 
For example, 



(A, 



-rS'\ 



{A)o = \^-rS'{^^ (4-5) 
and thus 

\^ = / rni -rs'\ ■ 4-^) 
(<T D e '^•^ )r 
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4b, r=0.25 I — ^ 
4b, r=0.5 ^--x- 

3d, r=0.25 -« 
3d, r=0.5 ^ B 



0.5 1 
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Figure 8: The index measured using the less important samphng for = 21. Set 4b should 
give 1 and 3d should give 0. 



The Wittcn index becomes 



w 



C 



(a[i^p]e-'-^p),,P 

^g(l-r)S^,-iE.M"</'?)^ 



for the configuration with the periodic boundary condition and 



(CT[i:'p]e(i-'')^APe- 



-s'- 



/r,AP 



(e(i-'')-5AP-5E.M='/'?)^^ 



(4.7) 



(4.8) 



AP 



is nothing but 



for the anti-periodic conditions. Note that in gauge theories changing 
changing the bare gauge coupling /3. 

Now the advantage of the less important sampling is clear. Let us set r = i for simplicity 
and consider the periodic case. In eq. (j2.10p the denominator has large fluctuations while 
the numerator has much smaller fluctuations. The behavior of the denominator requires 
large statistics and the unbalanced magnitude of the fluctuation with that of the numerator 
means a poor efficiency. In eq. (14.71) . the denominator is the square root of that of eq. (|2.10p . 
so the order of the fluctuation becomes much smaller. In addition, the numerator has the 
same order of the fluctuation. In total, we expect much more efficient measurement. 

In fig. |8l we show the result from the less important sampling method for periodic- 
boundary conditions. With a suitable choice of /z^ which gives a small error, we can identify 
the correct index as expected. Table [S] summarizes the result based on the less important 
sampling method with r = 0.5 and r = 0.25. To see the advantage, we plot the behavior 
of the error against the number of configurations in fig. [HI picking up set 4b. Contrary to 
the expectation, the magnitudes of the error are not so different from one another for larger 
statisticJ^. But the less important sampling case gives a more stable behavior. 

An interesting byproduct is that it makes it easier to fiip the sign ct[Z3p] for the n — S 
case. The less important sampling method reduces the hight of the potential barrier between 
the <t[Dp] — 1 region and the —1 region. This effect gives less autocorrelation for the n = 3 
case compared with the r = cases before (with either periodic or anti-periodic conditions) . 



^■^ If we had carefully optimized the value of /^^ , 
smaller than that for r = 0. 



the magnitude of the error for r ^ might be significantly 



15 



10^ 10^ 10" 10^ 10^ 

number of configurations 



Figure 9: Behavior of the error versus number of independent configurations for different 
value of r. They are for set 4b, A'' = 21 and /U^ with the smallest error. Dashed lines are 

10 X (TVconfs)-'/'. 



set 


N 


num. 




index 


set 


N 


num. 




index 


4b 


11 


99000 


1.0 


1.0049(21) 


4b 


11 


99000 


1.0 


0.995(3) 




15 


99000 


1.0 


1.007(8) 




15 


99000 


1.0 


1.010(11) 




21 


499000 


1.5 


1.018(16) 




21 


499000 


1.0 


0.987(17) 


3d 


11 


450 


3.0 


0.017(6) 


3d 


11 


12475 


2.5 


-0.003(4) 




15 


300 


3.0 


-0.008(6) 




15 


4975 


3.0 


0.005(4) 




21 


360 


2.5 


0.006(2) 




21 


9900 


2.0 


0.000(5) 



Table 6: Results from the less important sampling method with r = 0.25 (left) and r = 0.5 
(right), under the periodic boundary conditions, /x^ is chosen from 0.5, 1.0,. ..,3.0 to minimize 
the error. 
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5 Conclusion and Discussion 



In this paper, we proposed a method for measuring the Witten index using a lattice sim- 
ulation. The index is given by the average of the phase of the fermion determinant and 
we gave a suitable normalization to relate this average to the index. Our requirement for 
the lattice model is that it has at least one exact supersymmetry Q with — Q and the 
action is Q-exact. We checked that the method reproduces the known Witten index for a 
simple supersymmetric quantum mechanics. The application to higher dimensional models 
is trivial in principle. 

A disadvantage of the method is that it does not follow the philosophy of the important 
sampling. Therefore it requires high statistics. The situation would be worse the more 
degrees of freedom the model has. The efficiency for the higher dimensional model would 
be less practical. A less important sampling method may cure the situation to some extent. 

In addition, having more fields means having more parameters to be tuned to achieve 
better efficiency. The supersymmetric quantum mechanics we numerically treated has only 
1 parameter /x^, with our choice of Gaussian regularization functional. For the gauge theory, 
typically we have 3 such parameters; for scalar field, gauge field, and pseudo fermion field. 
One could even tune the regularization functional to have more parameters. 

For the gauge theory, however, it is difficult to give a suitable gauge invariant regular- 
ization functional of the gauge fields. We also have to normalize the gauge volume. In 
the following, we discuss a possible treatment on the basis of the gauge fixing, though it 
contains some subtleties. In the lattice simulation, we can fix the gauge in the measurement 
setting as many link variables as possible to unity using the gauge transformation. The 
normalization constant C should be obtained using only the remaining degrees of freedom. 
Note that in the constraint system the relevant degrees of freedom which contribute to the 
path integral are the remaining ones after the constraint is solved. Setting the link variables 
unity using the gauge transformation corresponds to solving the constraint. An important 
subtlety here, which may spoil this gauge fixing argument, is that it is not clear whether 
the gauge is completely fixed or not. Another subtlety is the regularization functional for 
which throughout in this paper we used a Gaussian functional. One can pick up the gauge 
field A" out of link variables and use the same Gaussian functional for , but this relation 
contains some lattice artifacts. Therefore, it should give the correct normalization only after 
taking continuum limit. In this procedure the gauge fix is crucial since otherwise non gauge 
invariant Gaussian functional cannot give a non-trivial expectation value. 

Our method can be presumably applied to a lattice model with a Q-closed term [16j . 
which has an exact Q-invariance. In this case we cannot exclude contributions from E > Q 
states to the index. We have to take Lphys oo limit for the temporal direction. 
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A Derivation of the path integral 

We give a short review of the derivation of the path integral. 
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A.l Bosons 

The bosonic Hamiltonian is 



H{p,q) = 



P 



+ V{q), 



2mQM 

where tuqu is a parameter with mass dimension 1 and V{q) is a potential. 
We calculate the amplitude 

{qN{t = -IT = -iT)\qo{t = -r = 0)) = {qN\e-^"\qo), 



(A.1) 



(A.2) 



where |go) and Iqn) are the initial and final state, respectively. The coordinate eigcn state 
\q) is normalized as 

{q\q')=5{q-q') 



and satisfies the completeness relation 



/oo 
dq\q){q\ 
-OO 



1. 



We divide the period T into A'^ small periods: 

T = aN. 

Inserting N — 1 completeness relations, we obtain 



(A.3) 



(A.4) 



(A.5) 



(9Jv|e \qo 



dqk]{qN\e \qN-i){qN-i\e ''"\qN-2) ■ ■ ■ {qi\e ""{qo), 



(A.6) 



and assuming a is small, we obtain 

{qk+i\e-''"\qk) = - aH)\qk) 

■ exp 



^ f dpk+i 
~J-oo 27r 



Pk+i ir ( qk+1 + qk \ , ■ I N 



exp 



J^QM / qk+\ - 



2 V 



• (A.7) 



Here, we have assumed H is Weyl ordered. To calculate (gfe+ij-ff j^fe), we have also used the 
momentum eigen state |p) which satisfies 



/OO 
dp\p){p\ 
-00 



1. 



(A.E 



roo /N-1 



In total, we obtain 

iV-] 

E 



X exp 



fe=0 



(A.9) 
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In particular, the partition function becomes 



Tr(e-^^) 



dqo{qN^qo\e \qo) 



N-l 

n 

fe=o 
X exp 



/ ™qm \ 
V 27ra / 



E 

fe=0 



dqk 

™QM / qk+1 — qk 



— oo 



V 



qk+1 + qk 



(A.10) 



where the periodic boundary condition q^ = qo is assumed in the exponent. Defining the 
1-dimensional field 4>k = y'WQM qk (and suitable rescalings for parameters in the potential) , 
we obtain 



N-l 



Tr{e-"^) = n 



k=0 



X exp 



— y r 



d(j)k 



N-l 
fe=0 



1 / <pk+l 



+ V 



4>k+l + <t>k 



(A.11) 



where the first line provides the normalized measure. The exponent goes to the (minus) 
action in a ^ limit. Note that the above derivation gives one specific lattice action at 
finite a, which may or may not be useful for the lattice simulation. 
For a dimensionless lattice field 



6|t^* = a ^<f>k, 



the measure becomes 



/^r=n[(i)7 

■' fc=o L ^ ^ " 



(A.12) 
(A.13) 



A. 2 Fermions 

We use the following convention for the Grassmann integration: 



/ 



dtp* dtp tptp* 



+1. 



(A.14) 



The creation and annihilation operators -tjj and -tjj^ satisfies the following anti-commutators: 



The Hilbert space is 2-dimensional and spanned by |0) and |1) as usual: 



^|0) =0, 

(0|0) = (1|1) = 1, 

The Hamiltonian for fermions is 



VS^IO) = |i), 

(0|1)=0. 



(A.15) 



(A.16) 
(A.17) 



H{'ii}\tp) = tp'^Mtp. 



(A.18) 
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It is convenient to introduce the following coherent state: 

^e-^*>) = |0)-V|l), 



where ip (without is a Grassmann odd c-number. This state satisfies 

m=m {=m)- 

The orthogonality and the completeness relations are 

(VI V') = l + V*V'' = e^*'^', 

dV*dV |V)e-'^*^(V| = |0)(0| + |1)(1| = 1. 

Inserting A'^ — 1 complete sets and using 

(Vfc+i|e-«^|Vfc) = (Vfc+il (1 - aH{rk+i,^k)) IVfc) 
= exp [-ai?(V'fe+i, Vfc) + V'fc+iV'fc] , 

we obtain 



(A.19) 
(A.20) 

(A.21) 
(A.22) 

(A.23) 



(Viv|e-^^|Vo) = n / '^^fc'^^'^ 



>fe=i 



N-l 



(A.24) 



fe=i 

+ V^Viv-i - aH{il)l,ipo) 



The trace of an operator A in the coherent representation is 

Tr(i) = - j d-ip* dip e'''''''{i)\A\il)) = (0|i|0) + (l|i|l). 



(A.25) 



Then we obtain 

Tr(e-^^) = - 



dr^ d,Poe'^o^'{Me-"^\^o) 



(A.26) 



= - / #fe#fejexp -^Xj |^^'(V'fe,V'fc-i) + ^V'fc(V'fc-V'fc-i)| 



- aH{tpQ,'ipN-i) +'>Po'^N-i + V'oV'o 



. (A.27) 



Replacing Vo — ^ — Vo ^'^d imposing the anti-periodic condition Vjv = — Vo and Viv = ~V'o! 
we obtain 



Tr(e-^^) 



/ dVfc dijjk exp 

; — n J / 



.k=0 



N-l 



fe=0 



a XI ^ ^ii^k+i^i^k) + ^tpl+iiipk+i - tpk) 



(A.28) 



In the continuum limit a ^ 0, the exponent gives the Euclidean action with anti-periodic 
boundary conditions. 

To conclude, the measure for the fermion field is 



N-l 

/ vr-Di>= n / drkdiPk 

•' k=0 •' 



(A.29) 
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